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A theory of traffic-measurement errors for loss systems with renewal 
input is developed. The results provide an accurate approximation for the 
variance of any differentiable function of one or more of the following basic 
traffic measurements taken during a given time interval: 
(i) The total number of attempts (peg count) 
(ii) The number of unsuccessful attempts (overflow count) 

(Hi) The usage based o?i discrete samples (TUR measurement) or on 
continuous scan. 

The approximation is given in terms of the individual variances and 
covariance functions of the three measurements. Asymptotic approxima- 
tions for these moments are obtained using the concept of a generalized re- 
newal process, and are shown to be sufficiently accurate for telephone 
traffic-engineering purposes. 

As an application of the theory, we examine the variances of the standard 
estimates of the load and peakedness (variance-to-mean ratio) of an input 
traffic stream for a time interval of one hour. Other possible applications to 
Bell System trunking problems are discussed. 

I. INTRODUCTION 

In the Bell System, there are a number of traffic measurements 
which can be made on any given trunk group. For a standard time 
interval (0, Q of one hour, the three most important measurements are : 

(i) A(t), the number of attempts (peg count); 
(ii) 0(t), the number of unsuccessful attempts (overflow count); 

and 
(Hi) L d (t), an estimate of usage based on 36 discrete samples (TUR 
measurement). 

When all three measurements are available, several statistics can be 
formed to estimate traffic parameters of interest. For instance, the 
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ratio 0(t)/A(t) is an estimate of call congestion. Two other important 
parameters are the peakedness (variance-to-mean ratio) and the load 
of the input traffic. An estimate of the load is given by the function 

_ Ld(Q/36 

A(t) 

while an estimate of peakedness is a complicated function of A(t), 
0(t), and L d (t) which is usually obtained by iteration using the Equiv- 
alent Random method. 1 

Since the trunk-engineering procedures are based on such estimates, 
it is important to know their statistical accuracy. For instance, it 
would be useful to know the error inherent in a prediction of the re- 
quired size of a trunk group (to obtain a specified grade of service) 
based on the estimates of offered load and peakedness of the input 
traffic. Such a result could be used to determine the number of single- 
hour measurements necessary to ensure a desired accuracy in the pre- 
diction, to determine the optimum number of measurements from a 
cost-effectiveness point of view, or to evaluate the consequences for 
trunk provisioning of using a given number of measurements. 

Many results concerning the accuracy of the individual traffic mea- 
surements (i) through (til) have been obtained previously, but most of 
these have assumed the arrival process to be Poisson. For example, 
assuming Poisson arrivals, the variance of the usage measurement 
L d (t) was obtained by Bene§, 2 and the variance of the measured call- 
congestion 0(t)/A(t) was given by Descloux. 3 More recently, the vari- 
ance of 0(t)/A(t) was obtained for arbitrary renewal input by Kuczura 
and Neal. 4 The variances of some nonstandard traffic counts were 
considered by Descloux, 5 and numerical results were obtained for the 
case of Poisson input. 

Using the concept of a multidimensional renewal process, we develop 
a general theory of errors which provides an estimate of the variance of 
any differentiable function of the measurements (i) through (Hi). Con- 
sequently, our results can be used to answer many questions similar to 
those mentioned above. The variances of the estimates of offered load 
and peakedness of the input traffic will be derived in Section IV as 
examples which illustrate our general theory. 

Section II contains the derivation of an approximation for the vari- 
ance of a function of the three traffic measurements. The approxima- 
tion is given in terms of the individual variances of (i), (ii), and (Hi) 
and the covariance functions between them. 
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Section III contains the mathematical model used to derive the vari- 
ances and covariances. Section V contains a summary and an outline 
of other possible applications. 

II. STANDARD ERRORS OF FUNCTIONS OF RANDOM VARIABLES 

For completeness, we present those results from the theory of stan- 
dard errors 6 which are required below. Let £1, £2, and £3 be random 
variables and g a real-valued function. Assume that £, has a mean 
0*i fiK£i> £2, £3) has finite mean and variance, and g is differentiable at 
the point (6 h 2 , 03)- Using a Taylor series expansion we have, to first 
order, 

rffci h, tz) - 9(ei, e*, 9 a ) m £ (fr - Oi) |g.. (1) 

We are assuming that the observation period will be sufficiently large 
so that the contribution of the higher-order terms can be neglected. 
This assumption will be justified in our model. 

Taking the expectation of (1), we see that the mean value of g is 
approximately g(6 h 2 , #3) since #[£,- — 0,] = 0. We also have 

Var &(&, fr, *,)] 

+ 2 tal c -«,« + 2 if 2 ||cov [f! , f ,]. (2) 

After setting fc = A(t)/t, £ 2 = 0(t)/t, and £3 = L d (t)/SQ the above re- 
lation becomes the starting point for our theory of traffic-measurement 
errors. It approximates the variance of any differentiable function of 
the measurements in terms of their variances and covariances. 

In the next section we derive expressions for the required moments. 
These, together with the first partial derivatives of g, approximately 
determine the variance of g. If the function g is too complicated to be 
differentiated analytically, differencing may be used to approximate the 
partial derivatives. An example of this procedure is given in Section IV 
where we discuss the variance of an estimate of the peakedness of a 
stream of offered traffic. 
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III. MATHEMATICAL MODEL 

Consider a system of c servers serving customers whose arrival 
epochs constitute a renewal process. We assume that the interarrival 
times are independent and identically distributed according to the dis- 
tribution F having mean 1/X, and that the service times are in- 
dependent and identically distributed according to a negative-exponen- 
tial distribution with unit mean. If all servers are occupied when a 
customer arrives, he leaves and has no further effect on the system. If 
an idle server is available when a customer arrives, service begins 
immediately. 

Let (0, t~\ denote a time interval of length t which commences at a 
stationary point for the arrival process.* (Such a point is often said to 
be chosen at random on the time axis.) Let A(t) be the number of 
arrivals and 0(t) the number of blocked attempts in (0, Q. Finally, let 

L(t) = [' C{u)du , (3) 

where C(u) is the number of busy servers at time u, be the total usage 
in (0, f\. Note that L(t)/t is a continuous-scan estimate of the carried 
load. 

As was pointed out in Section II, the individual first two moments 
of A(t), 0(t), and L d (t) and the corresponding three covariance func- 
tions are sufficient to obtain an estimate of the variance of any func- 
tion of these measurements. By numerical experimentation, we 
found that the covariance functions Cov£A(t)/t, L d (0/36] and 
Cov\_0(t)/t } Ld(i)/36] are, for our purposes, well approximated by 
Cov [_A(t)/t, L(t)/Q and Cov [0(0 A L(t)/Q, respectively. However, 
the variance of L(t) can be significantly smaller than the variance of 
L d {t) so that, in general, we must use Var \_L d (t)~] in our applications. 

In the next two sections, we derive the individual and joint moments 
of A(t), 0(t), and L(t). In Section 3.3, we obtain the variance of L d (t). 

3.1 A Multidimensional Renewal Process 

Assume that the system described above is in statistical equilibrium*, 
let t n , n = 0, 1, 2, • • • , be the instant of time at which the nth overflow 



'That is, the time until the first arrival after t = has the remaining life-time 
distribution H(t) = X Jo' [1 - F(x)^dx. 

f That is, the system has been in operation sufficiently long prior to t = so that 
system-state probability distribution at t = is the limiting (or stationary) dis- 
tribution P|C(0) = k\ = p k = limu.,. P\C(u) = A:}. It follows that for any t ^ 0, 
P{C(t) = k\ = Pk, i.e., the process \C(t), i^Oj is stationary. 
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occurs, to < < ti < to < ■ ■ ■ , and define X„ = t„ — t n -i- Now let 
K„, n = 1, 2, ■ ■ ■, be the number of arrivals in (t n -i, t*] and 



/„ = f tn C(u)du, n = 1,2, 

J <n-l 



be the total usage in (t n -i, t„2- 

Since holding times are exponential and arrival epochs constitute a 
renewal process, the sequence of times t„, n = 0, 1, 2, • • • , are regenera- 
tion or renewal points in our model. Hence X», K„, and /„, n = 1, 
2, • • • , are sequences of independent and identically distributed 
random variables. 

If we now define the row vector 

x„ = (1, K,„ /„), n = 1,2, ■•■, 

then it follows that jx„, X„\, n = 1, 2, • • • , is a multidimensional 
renewal process. 7 Moreover, setting 

r,(t) = SX„ , 

where the sum is taken over all n such that < t n ^ t, we see that for 
large t, 

v (t) - (0(t), A(t), L{t)). 

Since this formulation corresponds to the concept of a generalized 
renewal process as communicated by J. M. Hammersley in the dis- 
cussion of W. L. Smith's paper, 7 his results apply directly to our model. 
In particular, we shall use his equations (25) and (26) to compute the 
moments of -q(t). 

Let fi„(c) = E{X1) be the nth moment of the interoverflow times 
from a group of c servers and 



v n = |" f»dF(£) 



For brevity, we denote the arrival intensity vi x by X. From Ref. 4, we 
have the first moments of A(t) and 0(t) already computed: 

ELA{t)~\ = Xt, 

t (4) 

From eq. (25) of Ref. 7, we have 
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where 

a> n (c) = #[/?]■ 

Again, from Ref. 4 and Ref. 7, the variances and covariances of the 

three measurements for large t, omitting terms which behave as o(t), 
are given by 

Var[A(*)]~4[> 2 -,n, 

Var[0«)]~-3^[M2(c)-/x?(c)], 

Var[L«)] - jj^j |rf(c)«.(c) + m 2 (c)co?(c) 

-2 Ml (c)a, 1 (c)£[Z 1 / 1 ]}, (6) 
Cov[A(0, 0(t)~\ ~ -i, (\ M2 (c) - ^[iTxZx]} , 

Gov [O(0, L(0] ~ to Imi(c)«i(c) - m(.c)E[XJ{\h 
Cov[A(0, L(0] ~ jj^ {«(e)*LWi] + MeWe) 

We now need to compute the various moments and joint moments of 
X\, K\, and I\ appearing on the right-hand side of (4), (5), and (6) 
in order to evaluate the approximate expressions for the moments and 
joint moments of A(t), 0(t), and L(t). Note that the mean and variance 
of A{t) are known since X = vf 1 and v* are computed directly from F. 

3.2 The Joint Distribution of K\, X h and Ii 

The development here parallels that of Section 2.2 in Ref. 4. Let 
h c (w, r, n) be the joint density function defined by 

d 2 
h e {w, r, n) = j^ P{X 1 ^w,I 1 ^r,K 1 = n}- 

By considering the two mutually exclusive events {the cth trunk re- 
mains busy throughout (0, w)\ and {its complement}, and using a 
renewal-type argument, we arrive at the following integral equation: 

n— 1 r w r u r r—v 

h c (w, r, n) = e~ w h c -i(w, r — w,n) + £ / / / e~ v h c -i(u, s, k) 

t=i J o J o Jo 

X h e (w — u, r — s — v, n — k)dsdvdu , (7) 
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in which the time variables u and v run concurrently from an overflow 
epoch. 

The following boundary conditions hold : 

h c (w, r, n) = 0, for r > cw or r < , 
h e (w, cw, n) = e~ ew f(w)6i, n , 

where f = F' and 8i, n = 1 for n = 1 and is zero otherwise.* 
If we define 

y c (x,y,z) = £ / / e- xw -» T h e (w, r, n)z n drdw , (8) 

n=l J JO 

then it follows from (7) that 

y ( X v Z ) = (y + i)7^-i(g + y + i, y, «) . , 9) 

7=^, y , Z ) y + l _ 7c _ i(x> yt z) + 7c _ i(:c + y + lf y> f) w 

Motivated by the work of Riordan 8 and the success of the approach 
taken in Ref. 4, we set 

where Do(x, y, z) = 1, and, as can be seen by setting c = in (8), 

where 

0(a) = /" " e- zt dF{t) • 
Furthermore, for m ^ 1, 
D m+1 (x, y, z) 

= D m {x, y, z) + [|^y - 1 ] D m (x + y + 1, y, g) • (12) 
If we now define 

x '- x ^** ) - 1 - «(»+** + a ' (13) 

then using (11) and (12) and mathematical induction one can show that 
D m (x,y,z) = 1 + f(- l)'(*)x*i---X*-r (14) 

* In our model we assume that the interarrival-time probability distribution func- 
tion is differentiable. However, with more formalism, the same results can be obtained 
for the more general case; e.g., the one-point distribution function for constant 
interarrival times. 
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Now, since 






x = p = 
z = l 



for A; = 0, 1 and i, j ^ 0, we can compute the required moments di- 
rectly by means of differentiation. Omitting all of the details of the 
operations indicated, we obtain 

Wl (c) = D e - 1. 
where 

D c = D c (l,0, 1) = 1 + t (-D^Q^AoAx-.-A^x, 
with 

A* = 1 - 



4>{k + 1) 



Note that D c is the reciprocal of the generalized Erlang-B blocking 
probability B c . Moreover, with the aid of the results obtained in Ref. 
4, we have 

«,(c) = 2[«i(c) + 1] ± l>0) + 1] - 2|7>< W0 > + Dr o) i 



y=i 



V2 

where 



mi(c) = - /«i(c) + 2 M i(c) £ miW - 2v 1 Z>i 100 >, (16) 



Ai+i+k 

1=2=1 



«~ = =S=B3 fl <<*. ». •) 



dx i dy t dz 

The required derivatives are given by 

«-.-j(-iy(;)w.A H [g + g+...+fe], 
w--£(-«'(;)i-.-i«g+S+- + fe]. 

where £2* is the derivative of \k(x, 0, 1) evaluated at x = 1, i.e., 

_, _ »'(fc + 1) 
iZt ~ 0«(* + 1) ' 

and fl* is the derivative of \*(1, y, 1) evaluated at y = 0, i.e., 

A#'(fc -f 1) 1 



Q* = 



2 (fc + 1) *(fc + 1) 
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Similarly, we have the joint moments 

ElXthl = 2*i[«x(c) + 1] t [o) 1 (j) + 1] 

- (1 + rJD?™ - Vl D^°\ 
E\_K x X x -\ = m,(c) + - mi(c) + t w(fc) + »iDf» - D< 100 \ (18) 

E\K x l{\ = 2[a, 1 (c) + 1] f [«i(c) + 1] 

_£(ioo> _ ^toio) + D m» i 
where 

D », ) = t ( _ 1)J (^ AoAi ... A ,_ i [i. + i i+ ... + _L i _,]. 

3.3 Variance of Discrete-Scan Estimate of Usage 

Our present mathematical model assumes that the measurement of 
usage, L{t), is made by means of continuous scanning, as can be seen 
from the definition of L(t) in eq. (3). In practice, however, usage is 
estimated by discrete scanning. The number of busy trunks is sampled 
at constant intervals of time, say r, and the integral in (3) is replaced 
by the finite sum 

Lait) = L Cikr), (19) 

k=\ 

where 

nit) = max ik:kr ^ t). 

This procedure introduces a sampling error in the evaluation of the 
integral. As we shall see later, the difference between the variances of 
Lit) and Ld(t) can sometimes be large enough that the discrete-scan 
variance of usage must be used to estimate accurately the variance of 
g in eq. (2). In this section we indicate how the variance of L d it) is 
computed. 

Let n = nit) be the number of discrete samples in (0, t~\. In trunking 
applications, t is usually taken to be one hour (about 20 mean holding 
times) and, since t is normally set at 100 seconds, n = 36. Since the 
process \Cit), t ^ 0} is stationary, from (19) we have 

Var [L d (0] = nR(0) + 2 £ (n - j)RiJr), (20) 

where R is the covariance function defined by 

Rit) = £[C(0)C(0] - £[C(0)]£[C(<)]. (21) 
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Since P{C(0) = k\ = P[C(t) = k} = p k , 

R(t) = £kp k ± jP kj (t) - ml (22) 

k-o y-o 

where 

P kj (t) =P{Cit) = j\C(0) = k], (23) 

and 

m, = £ kp k . (24) 

*— i 

The problem of determining the transition function in (23) has been 
treated by Takacs in Chapter 4 of Ref. 9. However, he uses a renewal 
point for t = 0, i.e., his origin is chosen at a point immediately after an 
arrival has occurred. His result, though not directly applicable, can be 
modified in a straightforward manner to take account of our different 
location of the origin. We state here the analogue of his Theorem 3 for 
the case of a stationary origin and give a proof in the appendix. We 
use 1/m to denote the mean service time throughout the statement and 
proof of the theorem. 

Theorem: Let t = be a stationary point for the arrival process described 
above. Then the Laplace transform of {28) is given by 

ir kj (s) -f* e-> l P kj {t)dt = £(-1)'-' Q) Me), (25) 

where 

1-0( 8 + im) **.■(«) , /fc\ 1 [, 0(« + i/O l 
PkiW ~ <f>(s + v.) (s + t/0 + \ i) (s + *m) L *(• + W J ' 

X U S ( i) CM \l S \j)c~M4>(s + jm) J 

c -»-4r^fft5. , = 0,1,2,.., 

C-x(s) = 1, 
</>(s) =j"e-*dF(f), 
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and 0(s) is the Laplace-Stieltjes transform of the distribution {24), that 
is, 

o 

Taking the Laplace transform of (22), substituting for wkj(s), and 
simplifying, we obtain the following expression for p(s), the Laplace 
transform of the covariance function R(i) : 

.. m 2 |\ a f 1 - <f>(s + 1) 1\ m? 

p(«) - jqrr J 1 - r+T [ ^ + i) J) - t 

+[^ + «£(;)^]{t^(;)a 

, Xm 2 Cp(a) «, /c\ 1 
+ (s + 1) d(«) ft VjV CX«) 



h (s + j)Ci(s) I <f>(j) ' 



-+-' -^J^^]}, (26) 



1 ~ »(J + 1) 

4>U + 1) 



where 



Co = i, 

and mi and m 2 , the first and second moments of the distribution {pk\, 
are given by 9 

W»l - E fcp* = - (1 - Be), 
A-i M 

c 

m 2 = T, fc 2 ?* = ™i + X[J5i — cB e ]. 

Note that Wi is the carried load and 2? e is the generalized Erlang-B 
blocking 

«•-/£(;) 4- 

Equation (26) has been inverted analytically for the case of Poisson 
input. 2 When </>(s) does not have the simple expression of this special 
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case, analytical inversion appears to be complicated. However, for the 
purpose of computing the variance of L d (t) for our trunking applica- 
tion, it is unnecessary to obtain an explicit inverse of p(s). We found 
that the numerical inversion scheme described by Eisenberg 10 is com- 
putationally efficient and gives satisfactory results. 

To illustrate the difference between continuous-scan and discrete- 
scan measurements for the case when the input traffic is of the overflow 
type, we computed the estimates of the variances of L(t)/t, the con- 
tinuous-scan estimate of the carried load, and L d (t)/36, the discrete- 
scan estimate, for various trunk-group sizes and t = 20 mean holding 
times-i.e., about one hour. For these results, the interarrival-time dis- 
tribution of the arriving traffic was obtained by using the Interrupted 
Poisson process with a three-moment match. 11 

The case for ten trunks is typic al and is pres ented in Fig. 1 where 
a Ld = War \_L d (t) /36] and a L = War \_L{t)/f\ vs a are graphed for 
z = 1,2, and 4. 

Since the variance of LdOO/36 must be at least as large as the vari- 
ance of L(t)/t, our results show that the asymptotic estimate of 
Var \_L(t)/f\ has a small positive bias, especially for low loads. Our 
simulation results verify this observation and also indicate that the 
asymptotic approximation becomes more accurate as the input load 
increases. Notice that the variance of L d (i)/3Q is about equal to the 
variance of L(t)/t at low loads. As the load increases, the relative error 
introduced by discrete scanning can increase substantially. Finally, we 
found that for fixed load and peakedness, the relative difference be- 
tween Var [Ld(0/36] and Var [L(t)/t] decreases as the trunk-group 
size increases (an effect not shown in the figure). 

IV. TWO APPLICATIONS 

We give two applications of our results, in which we obtain the ac- 
curacy of the estimates of two traffic parameters. In the first example, 
the parameter is the offered load as given in Section I. For the second 
example, we discuss an estimate of the peakedness of the offered traffic. 

4.1 Accuracy of an Estimate of Offered Load 

Suppose we have observations A(t), 0(t), and L d (t) recorded. Then 
for the measurement period (0, t], St, an estimate of the offered load (in 
erlangs), is given by 

& = g£A(t), 0(0, Um = ^Jffi ■ (27) 

l ~Mt) 
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Fig. 1 — Standard deviation of carried load measurements vs offered load using 
discrete-scan and continuous-scan measurements on a 10-trunk group for t = 20 
mean holding times. 

Obtaining the required derivatives of g as indicated in eq. (2), sub- 
stituting into (2), and simplifying, we have the following expression 
for the variance of the off ered-load estimate in (27) : 



Var [<$] = 



{ 



Var [0(0] + Var 



[^)](4) ! 



Hi - B c y \ 

+ B\ Var [A(0] + 2 Cov [£(*), O(0] 

- 2B C Cov [L(0, A(t)2 - 2B C Cov [O(0, A(t)] 



(28) 



Now using eqs. (6), (15), (16), and (18) to substitute for the various 
quantities on the right-hand side of (28), we can compute Var [a]. 

To test the approximation (28), we computed the variance of &, as 
outlined above, for trunk-group sizes of c = 10 and c = 40 trunks, for 
input traffic streams of the overflow type having different combinations 
of load and peakedness values. We also used a computer simulation to 
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Fig. 
trunks. 



2 — Standard deviation of offered-load estimate vs offered load for c = 10 
The measurement interval is 20 mean holding times. 



estimate Var [<$] at several points. The numerical results are displayed 
in Fig. 2 for c = 10 trunks and Fig. 3 for c = 40 trunks where 
<x a = War [VS] vs a is displayed for 2 = 1,2, 4, and 10 (again for t = 20 
mean holding times). 
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Fig. 3 — Standard deviation of offered-load estimate vs offered load for c = 40 
trunks. The measurement period is 20 mean holding times. 
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The simulation results indicate that for c = 10 and c = 40 the 
asymptotic approximation for Var £<$'] is quite accurate for all ranges 
of load and peakedness of interest in trunking applications. 

We obtained almost identical results for Var [<2] regardless of 
whether we used Var [L d (<)/36] or Var [L{t)/t~\. Apparently, for low 
loads (low blocking probability) the accuracy of & is dominated by the 
accuracy of the usage measurements while at high loads (blocking 
near 1) the accuracy of the call-congestion estimate is the dominant 
factor. Since the relative difference between Var [L d (<)/36] and 
Var [_L(t)/t] is small for low blocking probabilities, we see that the 
accuracy of d is not significantly affected by the TUR sampling error. 

4.2 Accuracy of an Estimate of Traffic Peakedness 

When all three of the measurements A(t), 0{t), and Ld(t) are avail- 
able for a final trunk group of c trunks in an alternate-route network, 
the peakedness z of the input traffic is estimated in the following man- 
ner: First an estimate of offered load 6. is determined as described in 
the preceding section. Then an estimate z of z is obtained by iterative 
methods (using the Equivalent Random method 1 ), such that a single 
overflow stream having load & and peakedness z would experience the 
call congestion 0(t)/A(t) or, equivalently, the resulting carried load 
would be L d (0/36. 

Thus, there is a well-defined procedure for determining a unique 
value for z corresponding to A(t) ^ 0(t) > and L d (t) > 0, i.e., we 
have the estimate 

[,4(0 0(t) L d (t)l 

*- g [—' — '-w\ 

in the required form. However, there is no explicit analytical expres- 
sion for g which can be used to obtain the derivatives needed in (2) to 
obtain the variance of 2. 

In such cases, it is natural to estimate the partial derivatives by first 
differences. For example, 



dXi 



glOdl + A), fl 2 , 3 ] - g[0 u 2 , jj 



where A is a small positive number. Numerical experimentation indi- 
cated that A = 0.001 gives sufficient accuracy for the present applica- 
tion. Using the first-difference approximations as illustrated in (29) 
for the derivatives in (2) we have an estimate for the variance of z. 

We computed the resulting approximation for c = 10 and c = 40 
trunks for a range of offered loads a, several values of peakedness z, 
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Fig. 4 — Standard deviation of peakedness estimate vs offered load for c = 10 
trunks. The measurement interval is 20 mean holding times. 

and t = 20 mean holding times. We also compared our approximation 
with results obtained by simulation. The resul ts are d isplayed in Figs. 
4 (for c = 10) and 5 (for c = 40) where a, = War [i] vs a is given for 
z = 1, 2, 4, and 10. Note that <r z /z -0.2 for large a, independent of c. 



ANALYTICAL APPROXIMATION 



X POINTS OBTAINED BY SIMULATION 




OFFERED LOAD.cv, IN FRLANGS 

Fig. 5— Standard deviation of peakedness estimate vs offered load for c = 40 
trunks. The measurement period is 20 mean holding times. 
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Fig. 6 — Comparison of the standard deviation of peakedness estimates using 
discrete-scan and continuous-scan measurements of carried load on a 40-trunk group. 



In general, the simulation results are in good agreement with the 
approximation. The curves are plotted either for a ^ z, or else for call- 
congestion exceeding 0.01, the range of interest for trunking applica- 
tions. When a is smaller than z or the call-congestion is much smaller 
than 0.01 (not shown in the figures), the value obtained from the ap- 
proximation for Var \z~\ occasionally tends to be larger than that ob- 
tained hy simulation. Hence, the approximation may not be adequate 
for such applications. 

In the preceding section, we noted that essentially the same results 
were obtained for Var \Jk\ regardless of whether we used Var [L{t)/t~] 
or Var [L d (£)/36] in the computations. In contrast, the variance of z 
is very sensitive to the variance of the usage measurements. That is, 
Var [L d (0 /36] is required to obtain an accurate approximation for 
Var [z]. The errors which can result from using Var \_L(t)/f\ instead 
of Var [L d (0/36] are illustrated in Fig. 6 for the case of c = 40 trunks. 
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V. SUMMARY AND OTHER APPLICATIONS 

5.1 Summary 

We derived an approximation for the variance of any differentiate 
function of the three basic traffic measurements-namely, peg count, 
overflow count, and usage (TUB.). The approximation is expressed in 
terms of the first partial derivatives or first differences of the function, 
and the individual variances and covariances of the measurements. 
Except for the variance of the TUR measurements, asymptotic ap- 
proximations for the required moments were obtained by an applica- 
tion of Hammersley's generalized renewal theory. 

The variance of the TUR was given as a sum involving the covari- 
ance function for the number of busy servers at equally spaced scan 
intervals. The Laplace transform of the covariance function was de- 
rived and inverted numerically using an inversion technique described 
by Eisenberg. 

The results were then applied to obtain approximations for the 
variances of estimates of offered load and peakedness (variance-to- 
mean ratio) of a stream of traffic of the overflow type submitted to a 
loss system. The approximations were in good agreement with results 
obtained by simulation. 

5 . 2 Other A pplications 

5.2. 1 Offered Load Estimates Based on Usage Measurements 

At present, A(t) and 0(t) are not always measured on primary high- 
usage trunk groups. Estimates of the single-hour offered load and call- 
congestion for such groups are obtained from the TUR measurement 
L d (t) with an iterative procedure based on the Erlang-B theory. Using 
the techniques presented above, one could compare the accuracy of 
these estimates with that which would be obtained by using all three 
measurements. It should then be possible to evaluate the difference in 
statistical accuracy that results from using (or not using) the additional 
measurements. 

5.2.2 The Optimum Number of Single-Hour Measurements 

Normally, 20 single-hour measurements are used to obtain estimates 
§ of the correct number of trunks s required to obtain a specified grade 
of service. For example, on final trunk groups, the 20 single-hour esti- 
mates of call-congestion 0(t)/A(t) and usage L d {t) are first averaged 
and then used to estimate an average load and peakedness of the input 
traffic. These average values are used to obtain L 
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It has been proposed that the number of measurements be reduced 
in order to lower the data handling costs. However, reducing the num- 
ber of measurements would increase the variance of §, i.e., decrease the 
accuracy of the trunk estimates. It appears that an optimum number of 
measurements could be determined by minimizing a function of the 

form 

e(N) = JC Var [*(#)] + C N , (30) 

where Cn is the cost of taking and processing N measurements, §(N) 
is the estimate of s based on N measurements, and 3C is a cost associ- 
ated with inaccurate trunk estimates. The precise form of the function 
might require modification. However, the basic idea is to trade off an 
increase in the accuracy of the provisioning process due to more ac- 
curate trunk estimates (as N increases) against a corresponding in- 
crease in cost. 

It appears that one can obtain an approximation for Var[£(iV)] 
using an extension of the ideas presented in Sections II and III in order 
to account for the effects of day-to-day variation in the offered loads. 
However, a realistic model to justify (30) or to obtain 3C and Cn will 
require further study. 

APPENDIX 

We prove here the theorem stated in Section 3.3. We shall need the 
following lemma. 

Lemma: For the model described in the text, let Y(t) be the number of busy 
servers at time t, t = be a stationary point, and Y(0) = i. Now let Y n 
be the number of busy servers found by the nth arrival and t n be the time of 
the nth arrival. For n = 1,2, • • • , and r = 0,1, • • • , c, define 



and 



Then we have 



AP(s) = E je-'« ( Y r n )| 7(0) = i\ , (31) 

*.>(s) - t A&. (32) 



' U £■ ( j) cMll S\;7^55*(« + jm)J 

"" L £ \j) cM \[ &1 \j) c~M *(« + jl) Jj ' (33) 
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where 

C-l(s) - 1, 



0(«) = f 
J 



e-* l dF(t), 



and 



0(s) = / " e-"dF(t) = x/"°V"[l - F(0]<« 



-£[1 -*<•)]■ 

o 

The proof below is essentially the same as the proof of Lemma 1 of Ref . 
9, modified to account for the stationary origin. For n = 1, 2, • • • , we 
have 

E { e -"» + l ( 7 ; +1 )|l r n = h tn + l ~U = x\ 

e -(.+r») X (i +1 \ E{e- St "\ Y n = j\, for j < c, 
fl -(H-r,). ^ c \ £{ e -.<«| r„ = c} for j = c, 

because under the given conditions F„ has a binomial distribution 
with parameters j + 1 (for j = 0, 1, • • • , c — 1) or c (for j = c) and 
e - "*. If we remove the condition t n+ i — t n = x, that is, multiply by 
dP{t n +i — t n ^ x\ and integrate over all x, we obtain 

tf L-.|.«( F ; +1 ) | y. = jj = *(« + r M ) ( j r + 1 )^{«—-|F. = j}, 

for j = 0, 1, • • • , c — 1 and 

EJe-*^( r j; +1 )|y n = c\ = </>(« + r M ) (^£{e-"»|F„ = c}. 

If we multiply the corresponding equations by P{Y n = j] and add 
them for j = 0, 1, • • • , c, then we get 

A& +1 \s) = *(« + r/0 [a&(s) + AftUOO " ( r 1 j ) *#«] . ( 34 ) 
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for r = 1, 2, • • • , c and 

A®(s) = 0(s)[0( S )]"-i. 

Since F(0) = i and t = is a stationary point, we have 

A\?(s) - (*)*(• + */•). 

Forming the sum (32) we get 



or 



M«) = 



4>(s + r/i) f / i \ 0(s -(- rju) 



o 



[1 - <t>(s + rn)ll\r / <f>(s + rp) 

+ h.r-i(s) -( r I x )*.(■)]■ (35) 
Dividing both sides of this equation by C r (s) we get 



tir(s) _ \fr,>-i(s) 



/A^s + tm) / c \ 



C r (s) Cr-l(s) C r -l(s) 

Adding these equations over r, r — 1, r — 2, ■ • ■ , 1 we obtain 

*■>(«) _ f /t W(a + j/i) 1 .. M y/ c \ 1 

C r («) ft \ J ) *(■ + jm) Ci-i(s) ^ ,cl ; ft V J - 1 / C-iOO ' 

t = 1, 2, • • • , c. (36) 
Setting r = c in (36) we get 

^ (s) " £ ( J ) WTJ5 -BOSS/ S, ( i) zfe- (37) 

Substituting (37) into (36) we obtain ^,>(s) for r = 1, 2, • • • , c. If 
r = 0, then 

M«) = £ ^WIX*)]"- 1 

n-1 
0(S) 



1 - 0(a) 
This completes the proof of the lemma. 
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We now prove the theorem stated in Section 3.3. Again, the proof is 
a slight modification of the proof of Theorem 3 in Chapter 4 of Ref. 6. 
Let us define the binomial moments 

BM = £ (*)*««), i, r = 0,1, ■■■ ,c. 

From the definition it follows that 

PM = ^ (-1)^(1) BM, (38) 

so that setting 

j8 t >(s) = r e-«B ir (t)dt 

J o 

and forming the Laplace transform of (38), we get eq. (25) of the 
theorem. It remains to determine /3,>(s). 

Let Y(t) be the number of busy servers at time t and let t — be 
a stationary point. The times between those successive arrivals which 
find j servers busy are independent and identically distributed ran- 
dom variables. Hence, the sequence of epochs immediately preceding 
those arrivals which find j servers busy constitutes a renewal process. 
If 7(0) = i, we will denote the renewal function of such an imbedded 
renewal process by Mij(t). 

It may be helpful to recall that Ma(t) is the expected number of 
those calls which arrive in the time interval (0, t~\ and find exactly j 
servers busy, given that initially there are i servers busy. Hence, we 
can write 

Mi&t) = t P{tn ^t,Y n = j\ 7(0) - i\, (39) 

n-l 

where t n and F„ have the same meaning as in the statement of the 
lemma. 

If no calls arrive in (0, t~\, then Y(t) has the binomial distribution 
with parameters i and e~"', and B ir (t) is given by 



(:) 



e -r„i[;i _ F(*)j, 
where 

F(t) = x/'n - F(x)]dx. 
J o 

If one or more calls arrive in (0, f], let the last call's arrival epoch be 
u and at that instant let the number of busy servers be j. Now Y(t) has 
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the binomial distribution with parameters j + 1 (if j = 0, 1, • • • , 
c — 1) or c (if j = c) and e~" ( ' -u) . Thus, together we have 

+ C Z x ( J t * ) /„' «~ r " (l ~ ,,) £ 1 ~ F C ~ «) Wi(«) 

+ ( C ) J' e-«-"<'- tt >[l - /P(i - u)]dil/, e (u). (40) 
If we introduce the Laplace-Stieltjes transform 
M0 (s) =r e-'tdMM, 

J 

then from (40) we have 

1 - 0(s + r M ) f / i \ 1 - 0(s + r M ) 
Pi * 8) (s + r M ) LWl -*(« + »» 

+ X( i t 1 ) Mo(s) + (r) M,c(S) ]' (41) 
From (39) we have 

mo(s) = t P\Y n = j)tf{«-'-|F. = h Y(0) = i\ 

n-l 

and hence by (31) and (32) we get 

- *•>(«). 

Thus, @ir(s) can be written in the following form 
s + rn l\r / 1 - 0(s + r/t) 



+ lfc.r-l(*) -( f * x)^*)] 



If we take relation (35) into consideration, then this formula can be 
simplified to 

1 - 0(8 + I7t) fr,(a) /i\ 1 [ 0(8 + r M ) I 

p<rW 0(s + r M ) (s + rn) ^ \ r ) (s + r/x) |_ *(* + r M ) J ' 

This completes the proof of the theorem. 
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